Physics of collisionless reconnection in a stressed X-point collapse 
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Recently, magnetic reconnection during collisionless, stressed, X-point collapse was studied using 
kinetic, 2.5D, fully electromagnetic, relativistic Particle-in-Cell numerical code [D. Tsiklauri and 
T. Haruki, Phys. Plasmas 14, 112905 (2007)]. Here we finalise the investigation of this topic 
by addressing key outstanding physical questions: (i) which term in the generalised Ohm's law is 
responsible for the generation of the reconnection electric field? (ii) how does the time evolution of 
the reconnected flux vary with the ion-electron mass ratio? (iii) what is the exact energy budget of 
the reconnection process, i.e. in which proportion initial (mostly magnetic) energy is converted into 
other forms of energy? (iv) are there any anisotropies in the velocity distribution of the accelerated 
particles? It has been established here that: (i) reconnection electric field is generated by the 
electron pressure tensor off-diagonal terms, resembling to the case of tearing unstable Harris current 
sheet studied by the GEM reconnection challenge; (ii) For mi/me 3> 1 the time evolution of the 
reconnected fiux is independent of ion-electron mass ratio; also, in the case of mi/me = 1 we show 
that reconnection proceeds slowly as the Hall term is zero; when mi/ me 3> 1 (i.e. the Hall term is 
non-zero) reconnection is fast and we conjecture that this is due to magnetic field being frozen into 
electron fluid, which moves significantly faster than ion fluid; (iii) within one Alfven time, somewhat 
less than half (~ 40%) of the initial total (roughly magnetic) energy is converted into the kinetic 
energy of electrons, and somewhat more than half (~ 60%) into kinetic energy of ions (similar to 
solar flare observations); (iv) in the strongly stressed X-point case, in about one Alfven time, a full 
isotropy in all three spatial directions of the velocity distribution is seen for super-thermal electrons 
(also commensurate to solar flare observations). 

PACS numbers: 52.35.Vd; 96.60.Iv; 52.65.Rr; 45.50.Dd; 96.60.pf; 96.60.qe 



I. INTRODUCTION 

In many astrophysical or laboratory plasma situations 
(a) plasma beta is small, indicative of large amounts of 
energy stored in a form of magnetic field and (b) there 
is a need to explain or provide plasma heating, as well 
as plasma particle acceleration. It is believed that in 
such situations magnetic reconnection, i.e. change of 
connectivity of magnetic field lines that penetrate the 
plasma, can serve as one of the important possible mech- 
anisms. There are different types of magnetic reconnec- 
tion. One of the key descriptors is plasma coUisionality, 
i.e. if plasma is coUisional then magnetic resistivity, t] (or 
more specifically 7]j term in the generalised Ohm's law) is 
responsible for breaking the frozen-in condition (enabling 
field line connectivity change). However, if plasma is col- 
lisionless, then other terms in the generalised Ohm's law 
may be more important. In this context, it is instruc- 
tive to look at typical spatial scales. Let us consider 
an example of solar coronal plasma. Typical width of a 
Sweet-Parker current sheet is given hy 5 = S~^^^L [l| (p. 
54). Fixing coronal temperature at 1 MK, Coulomb log- 
arithm at 18.0948, the Lundquist number (using Spitzer 
resistivity from Huba [1], p. 30) is 5.37933 x lO^^^ .^^j^ich 
for L = 10 Mm yields, S — 4.31157 m. Another way of 
looking at S is associating it also with the resistive length 
scale via Alfven time scale ta = L/Va, where Va is the 
Alfven speed (« 1 Mm s~^): i.e. S = S~^f'^L = r]TA/^J-o■ 
Typical scale associated with the Hall term in the gen- 
eralised Ohm's law at which deviation from electron-ion 



coupled dynamics is observed is, c/cupi = 7.20064 m (ion 
inertial length). Here particle density of n = 10^^ m^'^ is 
used. The fact that 5/{c/ujpi) < 1 points to a necessity of 
going beyond single fluid MHD approximation. The Hall 
term itself cannot break the frozen-in condition, its in- 
clusion into consideration ensures that the magnetic field 
is frozen into electron fluid. Below, we shall use this ar- 
gument to conjecture why reconnection is fast when the 
Hall term is included. Importance of different terms in 
the generalised Ohm's law is usually inferred by compar- 
ing the spatial scales associated with them to the resistive 
length scale 5. E.g. one of the other noteworthy scales 
is c/ujpe = 0.16804 m (electron inertial length) on which 
the electron inertia term operates. When 6/{c/ujpe) < 1 
then electron inertia would dominate over resistive diffu- 
sion [ij (p. 200). As can be seen from the above esti- 
mates S/{c/Wpe) — 25.6576 » 1 in the solar corona, thus 
electron inertia effects seem to be negligible. However, 
electron inertia in contrary to the Hall term (as well as 
the electron pressure tensor) can break the frozen-in con- 
dition and thus change the magnetic field connectivity. 
Because, of the fact that with increase of T (hot plas- 
mas), 6 gets progressively smaller thus effects other than 
the resistivity should be included. Indeed, collisionless 
(non-resistive) reconnection has recentlyattracted con- 
siderable attention (see Birn and Priest \3\ for a review) . 

One of the first studies of magnetic reconnection is 
stressed, X-point collapse [4;] (also see Chap. 7.1 in Priest 
and Forbes 0). The latter was using resistive MHD ap- 
proach. We recently revisited the problem in the regime 



of collisionless reconnection V\ . In Ref . V\ we studied the 
magnetic reconnection during collisionless, stressed, X- 
point collapse using kinetic, 2.5D, fully electromagnetic, 
relativistic Particle-in-Cell numerical code. We consid- 
ered two cases of weakly and strongly stressed X-point 
collapse. Where descriptors weakly and strongly refer 
to 20 % and 124 % unidirectional spatial compression 
of the X-point, respectively. Amongst other interesting 
outcomes, we established that within about one Alfven 
time, 2% and 20% of the initial magnetic energy can be 
converted into heat and accelerated particle energy in the 
cases of weak and strong stress, respectively. However, 
open questions remained: (i) which term in the gener- 
alised Ohm's law is responsible for the generation of the 
reconnection electric field? (ii) how does the time evolu- 
tion of the reconnected flux vary with the ion-electron 
mass ratio? (iii) what is the exact energy budget of 
the reconnection process, i.e. in which proportion initial 
(mostly magnetic) energy is converted into other forms 
of energy? (iv) are there any anisotropies in the veloc- 
ity distribution of the accelerated particles? Here we fi- 
nalise the study of magnetic, collisionless reconnection of 
a stressed X-point by providing answers to these ques- 
tions. 



II. THE MODEL 

The numerical code used here is 2.5D, relativistic, fully 
electromagnetic PIC code, with the initial conditions the 
same as in our previous work [y]. For completeness we 
re-iterate key points: Magnetic field configuration is an 
X-point without a guide-field 
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where Bq is magnetic field intensity at the distance L 
from the X-point (L is the global system scale) . a is the 
stress parameter, which prescribes the initial strength of 
magnetic pressure that collapses the system, due to lack 
of restoring force j5|. Using fioj — V x _B, a uniform 
current is imposed in the z direction. 
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Electrons and ions have uniform spatial, and Maxwellian 
velocity distributions throughout the system. For a = 1 , 
magnetic field geometry is completely symmetric, jz cur- 
rent zero (see Eqs. ([I])-©), and thus, such magnetic con- 
figuration is stable. For a > 1, stressed X-point starts 
collapse in the x direction because of the absence of a 
restoring force, causing time-transient magnetic recon- 
nection. The main parameters of the standard simu- 
lation model are as follows. The length of the system 
in two dimensions is L^ = Ly = 400A (this is exclud- 
ing so-called ghost cells), where A = 1 is the simula- 
tion grid size corresponding to electron Debye length. 



Ad = vte/ojpe — lA {vte IS clectron thermal velocity 
and ujpe is electron plasma frequency). The global re- 
connection scale is set L = 200A. The number density 
is fixed at no = 100 electron-ion pairs per cell. Hence 
the total number is 1.6 x 10^ pairs. The simulation time 
step is ujpeAt = 0.05. lon-to-electron mass ratio is fixed 
at rrii/me = 100 (which is varied in Fig. 2 (subsection 
III.B) only). The electron thermal velocity to speed of 
light ratio is Vte/c = 0.1. The electron and ion skin 
depths are c/wpe — lOA and c/ojpi = lOOA, respectively. 
The electron cyclotron frequency to plasma frequency ra- 
tio is Wce/'jJpe = 1-0 for magnetic field intensity, B — Bq. 
This ratio is close to unity in the solar corona, while it 
is much bigger than unity in the Earth magnetosphere. 
The electron and ion Larmor radii are vte/i^ce = lA and 
vti/uJci = lOA, where vu is the ion thermal velocity. Ini- 
tial temperatures of electrons and ions are initially set 
the same, T^ = T^. At the boundary {B — Bq at the 
distance L from the X-point), the plasma (3 — 0.02 and 
Alfven velocity, Vao/c = 0.1. Naturally these vary across 
the simulation box as the background magnetic field is a 
function of x and y. 

The boundary conditions on EM-fields are zero- 
gradient and also, tangential component of electric field 
was forced to zero, while normal component of magnetic 
field was kept constant, both at the boundary. This en- 
sures that there is no change in magnetic flux through 
the simulation box, i.e. the system is isolated. When 
colliding with boundaries particles are reflected. Thus 
our boundary conditions ensure there is no magnetic in- 
flux or mass transport across the boundaries. It has been 
also confirmed that the total energy in the system is con- 
served during the simulations to a good accuracy. 



III. RESULTS 



Before we address the outstanding questions, we refer 
reader to Tsiklauri and Haruki [y] for a more detailed 
description of dynamics of EM-fields, currents, and par- 
ticles. It is our intention to focus on the outstanding 
questions here. In brief, the previous results can be sum- 
marised as follows: when a > 1, the stressed X-point col- 
lapses in the x direction due to the absence of a restor- 
ing force, and hence time-transient magnetic reconnec- 
tion occurs. The fast reconnection regime is achieved. 
Initially uniform out-of-plane current becomes localised, 
peaking at many times its initial value in a time transient 
manner. Also, out-of-plane quadruple magnetic field is 
generated. High energy part of the electron distribution 
function exhibits a power-law behaviour. Sizable fraction 
of initial magnetic energy is converted into other forms 
of energy. 
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FIG. 1: Line plots of different terms in the generalised 
Ohm's law along y direction, in x = 0, at ujpet — 170 
(time-transient reconnection peak) for a = 1.20. Solid 
lines in panels (a)-(d) show the different terms as follows: 
(a) VeyBx, (b) -{di/n)dPexz/dx, (c) -{di/n)dPEyz/dy and 
(d) —di{me/mi)veydvez/dy. The reconnection electric field 
Ez{0, y) (normalised to VaoBq) is shown with the dashed line 
in all four panels. Here y is normalised to c/ujpi, and thus 
varies between —2<y<2. 



A. Source of the reconnection electric field 

In order to understand details of the coUisionless re- 
connection process, we now focus on the question: which 
term in the generahzed Ohm's law is responsible for the 
generation of the reconnection electric field? We adopt 
an appoach used by Pritchett Q • The generalized Ohm's 
law can be written as (e.g. [1] p. 108) 

^ ^ ^ V • Pp nip ,dvp ,^ _^ ^ X 

E^-VeXB i_^^+ i/^.Vi/e , 3 

Uee e ot 

where E and B are electric and magnetic fields, v is 
plasma velocity, P is pressure tensor (3x3 matrix), n 
is plasma number density, m is mass and e is electric 
charge. The subscript e refers to an electron. Normal- 
ising space coordinates by global reconnection scale L, 
fluid velocity by Alfven speed Va, time by Alfven transit 
time ta{= L/Va), magnetic field by JSq, number den- 
sity by no and pressure tensor by Bq/ ^q, a dimensionless 
version of Eq. ^ can be obtained 

E = -vp X B - d,^-^ - d,^(^ + {vp ■ V)vp) , (4) 
Up rrii ot 

where di is the normalised ion skin depth (di = c/cUpiL). 
Note that strictly speaking we should have used tildes in 
Eq.(j4]) to denote dimensionless quantities, but we omit 
them for brevity. 

Let us focus on the out-of-plane component of the elec- 
tric field Ez at the magnetic null, which is a measure of 
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where d I dz 
model. 



is assumed because of 2D reconnection 



The pressure tensor is defined as Py — 
m J v^Vjf{r,v,t)dv, where m is mass, v is random 
velocity, the subscript i and j denote the components x, 
y or z, /is the particle velocity distribution function, 
f is position, and v is velocity. In order to get the 
pressure tensor, number density is calculated first, from 
n{r,t) = J f{f,v,t)dv. Mean velocity is also obtained 
via Vifjt) — (l/n) J vf{f,v,t)dv. For pressure tensor 
calculation, the number density, n and the mean velocity, 
V{f,t), is calculated by counting number of individual 
particles per cell and by computing the average velocity 
in each cell, respectively. We then estimate the random 
velocity, v' = v — V, which is used in the above definition 
of the pressure tensor Py. In PIC simulations, in 
practise, the summation of mv^V: over all individual 
particles is used. 

Figure [1] shows ^/-profiles of different terms in the gen- 
eralised Ohm's law at x = 0, for ujppt = 170 (time- 
transient reconnection peak). Here a = 1.20. Solid lines 
in panels (a)-(d) indicate different terms as follows: (a) 
VpyB^, (b) -{di/n)dPpj:z/dx, (c) -{di/n)dPpyz/dy and 
(d) —di{mp/mi)Vpydvpz/dy. The reconnection electric 
field Ez{0,y) is shown with the dashed line in all four 
panels. A boxcar average scheme with a width of 7 mesh 
points is applied for smoothing data. The other terms in 
Eq. ([5]) are negligibly small. Fig. [T] is analogous to figure 
5 from Pritchett [7] . Fig. [Ija) shows that in all regions 
except the magnetic null, (0,0), contribution to Ez{0,y) 
from the VpXB term is significant. However, VpyB^ is zero 
at the X-point (the magnetic null) . As seen in Fig. [H^b- 
c) , the off-diagonal components of the electron pressure 
tensor are major contributors to £'^(0,0). The electron 
inertia term also generates the electric field away from 
the X-point (see Fig.[T](d)). Thus, we conclude that the 
reconnection electric field is generated by the electron 
pressure tensor off-diagonal terms; and hence the latter 
are responsible for breaking the frozen-in condition. A 
similar conclusion was reached by Pritchett Q . This co- 
incidence seems somewhat unexpected, because X-point 
collapse considered here and onset of tearing instability 
considered by Pritchett Q are physically different. Sim- 
ilarity of the cause of breaking of the frozen-in condition 
in both cases can only point to a universal nature of this 
mechanism. 




FIG. 2: Time evolution of the magnetic flux difference be- 
tween the O and X points (i.e. amount of reconnected flux). 
The solid lines with progressively increasing thickness show 
cases of mi/m^ = 9, 25, 64 and 100, respectively. The dotted 
line shows electron-positron plasma case {mi/me = 1). The 
magnetic flux difference Atj) is normalised by (Boc/cupi) and 
then a unity is subtracted to start from zero. Time is nor- 
malised by the ion cyclotron frequency uid = eB/rrii. Here 
a = 1.2. 



B. Effect of variation of the ion-electron mass ratio 
and conjecture of fast reconnection 

The next question we consider is: how does the time 
evolution of the reconnected flux vary with the ion- 
electron mass ratio? Such question historically was rele- 
vant because of the inability of performing realistic ion- 
electron mass ratio (1836) numerical simulation, due 
to lack of computational resources. Although within 
our reach computationally, we do not show here results 
rrii/me > 100 because the total energy conservation er- 
ror (which is defined as {E{ujpet = 250) — E{ujpet = 
0))/E{ujpet = 0) and is indicative of the code accu- 
racy) starts to deteriorate to values of circa 10% for 
mi/me = 400, while for rrii/me = 100 it is 0.04% (both 
for Of = 1.2). In order to be able to compare our results 
with the previous work [8|, when varying mi/me, we ac- 
cordingly adjust number of spatial grid points and total 
time step. Such adjustments insure that spatial scale of 
the simulation box, is L^ = Ly = Ac/cupi, and the time 
scale, LOcit = 25. Thus when setting mi/me = 1,9,25,64 
and 100, Accordingly, the system size is adjusted to 
40, 120, 200, 320 and 400A. The global reconnection size 
is fixed at L = 200A. Ion cyclotron frequency for each 
case is defined using by the magnetic intensity at the 
boundary. 

In 2D the magnetic flux function can be defined as 
Tp = — J Bxdy = J Bydx. In our simulation, X-point is 
located at the centre of system {x^y) = (0,0), while O- 
points are at {x,y) = (0,-2) and (0,2). Note that spa- 
tial coordinates here are normalised by ion skin depth. 



and L.J., — Ly = Ac/ujpi. Therefore we can use the same 
definition of the reconnected flux as in the case tearing 
mode- unstable Harris current sheet Q. Figure [2] shows 
time evolution of the magnetic flux difference between 
O and X points (reconnected flux) for different ion-to- 
electron mass ratios, rrii/me = 1,9,25,64 and 100. We 
gather from this graph that time dynamics of the recon- 
nected flux does not depend on mi/rrie when rrii/me 3> 1 
and that reconnection is fast. In fact, the time derivative 
of the reconnected flux is the reconnection rate. Thus, 
the conclusion is that the reconnection rate is indepen- 
dent of the mass ratio (when mi/rrie 3> 1). As with above 
conclusion (in previous subsection) that the reconnection 
electric field is generated by the electron pressure ten- 
sor off-diagonal terms; again similarity with the tearing- 
unstable Harris current sheet holds, i.e. Hesse et al. [8| 
came to the same conclusion in their case. 



As a further test, we performed a numerical run with 
mi/rrif, = 1 (case of electron-positron plasma). One of 
the main conclusions of Birn et al. [3] was that as long 
as Hall term is included, the reconnection is fast. i.e. 
when electron and ion dynamics can be distinguished. 
They showed that slow reconnection occurs only in the 
case of single fluid resistive MHD (in which there is no 
distinction in the electron- ion dynamics). However, in 
two-fluid MHD or PIC simulation it is possible to switch 
off the Hall term by setting mi/rrie = 1 as this will make 
electron-ion dynamics indistinguishable. The result is 
given by the dotted line in Fig. ^ . It can be clearly seen 
that the amount of reconnected flux grows very slowly 
with time, indicating that the reconnection is slow, as 
expected. 



We propose the following conjecture to explain why 
the reconnection is fast when the Hall term is included. 
Inclusion of the latter means that in the reconnection in- 
flow magnetic field is frozen into electron fluid. As it was 
previously shown in Tsiklauri and Haruki (see their 
Figs. (7) and (11)) speed of electrons, during the recon- 
nection peak time, is at least 4-5 times greater than that 
of ions. This means that electrons can bring in / take out 
the magnetic field attached to them into / away from the 
diffusion region much faster than in the case of single fluid 
MHD which does not distinguish between electron-ion 
dynamics. In fact, in Fig. (2) the amount of reconnected 
flux attained by ojdt — 25 in the cases of rrii/me 3> 1 
and mi/rrie = 1 has the same ratio (0.11/0.03 w 4) as is 
the ratio of electron and ion speeds (« 4 — 5). 



It should be mentioned that although the importance 
of the Hall term for providing the fast reconnection has 
been firmly established, some results indicate [lOl. [ill. [l2| 
that the fast reconnection without it is still possible. 
This indicates that the issue of which physical factor(s) 
uniquely guarantee the fast rate is still open. 



(a) 




150 200 250 



50 100 150 200 250 




00 150 200 250 



FIG. 3: Time evolution of (a) magnetic field energy, (b) elec- 
tric field energy, relativistic (c) electron and (d) ion kinetic 
energies of the whole system for a = 1.20. These energies are 
normalised by the initial total energy. Time is normalised by 
the electron plasma frequency tjpe. 
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FIG. 4: The same as in Fig. El but for a = 2.24. 



C. Energy budget of the reconnection process 

The next question we address is: what is the exact 
energy budget of the reconnection process, i.e. in which 
proportion initial magnetic energy is converted into other 
forms of energy? 

Figure [3] shows time evolution of (a) magnetic field en- 
ergy, (b) electric field energy, relativistic (c) electron and 
(d) ion kinetic energies of the whole system for a = 1.20. 
According to the previous results [6], in this case the 
normalised reconnection rate peaks at i?^ = 0.11 at time 
Wpe^ = 170. Initially magnetic field energy is dominant, 
which constitutes 96% of the total energy of system. The 
rest 4% goes to the initial electron and ion kinetic ener- 
gies because we impose a non-zero current jz at i = 
according to Eq.® (as a > 1). We gather from Fig. [31^ a) 



that as the reconnection proceeds magnetic field energy is 
converted into other forms of energy. As it is also stated 
in Tsiklauri and Haruki [y] at lOpei- ~ 250, which corre- 
sponds to about 1.25 Alfven times, (0.96-0.94)/0.96 = 2 
% of the initial magnetic energy is released. Here we ex- 
plore partition into which other forms this energy goes 
into. Panels (b)-(d) in Fig. [3] show that all other forms 
of energy increase as time progresses. In particular, elec- 
tric field energy that starts from zero, attains value of 
0.0024, i.e. (0 - 0.0024)7(0.96 - 0.94) = 12% of the 
consumed magnetic energy. One can conjecture that 
ultimately this energy will go into particle kinetic en- 
ergy (as particles would be easily accelerated by electric 
fields). Relativistic kinetic energy of electrons attains 
(0.0205 - 0.0145)7(0.96 - 0.94) = 30%, while the same 
for ions (0.037 - 0.0255)7(0.96 - 0.94) = 58%. Given 
that electrons have small inertia and thus are more influ- 
enced by the electric field, we conjecture that within few 
Alfven times electron-ion kinetic energy partition (as the 
percentage of consumed magnetic energy) will be roughly 
40% - 60%. It should be mentioned that, in general, it is 
not easy for an electric field to accelerate particles, un- 
less it is parallel to the magnetic field or the magnetic 
field is small, such that the particles are non-adiabatic 
(loosely tied to magnetic field lines). The latter condition 
is more stringent for electrons. Movies of electric field in 
our numerical simulation show complicated, and yet co- 
herent, oscillatory patterns. We have not performed de- 
tailed analysis of identification of nature of these waves, 
but based on previous experience, the case without the 
guide field considered here, would support excitation of 
whistler waves as the X-point collapses. In turn, we con- 
jecture that the whistler waves are ultimately responsi- 
ble for the particle acceleration. More detailed analysis 
of this is needed in the future. Returning to the issue of 
the established 40% - 60% energy partition, Emslie et al. 
[131] showed that the energy of accelerated electrons is 
comparable to that of accelerated ions. However, they 
admit to large uncertainties in the ion energy spectrum. 
Despite of this, our simulation results broadly agree with 
the solar flare observations [13j . 

Previously we also considered strongly stressed X-point 
(a = 2.24) Q. In this case in 1.25 Alfven times, 
(0.9 - 0.72)70.9 = 20% of the initial magnetic energy 
is converted into other forms of energy (this is equivalent 
of (0.9 — 0.72) = 18% of the initial toiaZ energy; and as we 
saw in the weakly stressed case, the difference between 
the two is negligible. It is only with the increase of a 
the difference between initial magnetic energy and initial 
total energy becomes noticeable, because stronger initial 
currents (i.e. initial kinetic energy of particles) need to 
be imposed according to Eq.®). This 18% decrease in 
the magnetic energy is also corroborated in panel (a) in 
Fig. m Exact break down (partition) of the latter is as 
follows (based on panels (b)-(d)): electric field energy 
that starts from zero, peaks and then settles at 0.006, i.e. 
(0 - 0.006)7(0.9 - 0.72) = 3% of the consumed magnetic 
energy. Relativistic kinetic energy of electrons attains 
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FIG. 5: (a-c) Electron and (d-f) ion velocity distribution 
functions in x, y and z directions near the current sheet 
at the initial stage i = (dashed line), the peak recon- 
nection stage t — 170 (dotted line) and the final simula- 
tion time t = 250 (solid line) for a = 1.20. As in Ref. 6], 
here data are produced using the region of the current sheet 

{-2{c/ujpe) < X < 2{c/ujpe),-8{c/0Jpe) < V < 8(c/LJpe)). fe 

and fi are the number of electrons and ions, respectively. Ve- 
locity and time are normalised by light speed c and ujpe, re- 
spectively. 



(0.075 - 0.002)7(0.9 - 0.72) = 41%, while the same for 
ions (0.2 - 0.1)7(0.9 - 0.72) = 56%. As in the weakly 
stressed case, within 1.25 Alfven times, electron-ion ki- 
netic energy partition (as the percentage of total energy, 
which for a solar flare would be the total energy released 
by flare) is roughly 40% - 60%. This again is in accord 
to solar flare observations Emslie et al. [l3| . 



D. Properties of velocity distribution of the 
accelerated particles 

The final question we address is: are there any 
anisotropies in the velocity distribution of the acceler- 
ated particles? This question naturally comes to one's 
mind due to a recent study of Kontar and Brown [14| . 
who surprisingly found near-isotropic electron distribu- 
tions in solar flares, which contrast strongly with the ex- 
pectations from the standard model that invokes strong 
downward beaming, including the coUisional thick-target 
model. 

As in Tsiklauri and Haruki p|, here we consider two 
cases of the weakly and strongly stressed X-point. The 
results are shown in Figures O and [SI The following ob- 



FIG. 6: (a-c) Electron and (d-f) ion velocity distribution 
functions in x, y and z directions near the current sheet 
at the initial stage t — (dashed line), the peak recon- 
nection time i = 45 (dotted line) and the final simulation 
time t = 250 (solid line) for a — 2.24. Data are produced 
using the region of the current sheet {—l{c/ujpe) < x < 
Kyc/uipe), — 16(c7cjpe) < y < IQ^c/iOpe)). The normalisation 
is same as in Fig[S] 



servations can be made: 

In the weakly stressed case, for electrons we see appear- 
ance of super-thermal electrons towards the end of sim- 
ulation time (shortly after the peak of time-transient re- 
connection) mostly in y and z velocity distribution func- 
tion components. Dynamics of the flows and currents is 
presented in detail in Ref. _6J . Here we only mention that 
the reconnection inflow is in x direction, while the outflow 
is in y direction. Thus based on panels (a)-(c) in Fig. ([51) 
we gather that accelerated electrons (focus on solid and 
dotted curves) are due to reconnection outflow (in re- 
connection plane) as well as out-of-plane flow (which is 
triggered by the out-of-plane electric field generated at 
the magnetic null). For ions, at later stages of the recon- 
nection, in panels (d)-(f) in Fig. ([5]) we see (focus on solid 
and dotted curves) a superposition of two Maxwellian 
distributions in both reconnection inflow (along x) and 
outflow (along y). These seem to be created by reconnec- 
tion flow dynamics. In z-direction we see a shifted (also 
somewhat broadened by the heating) Maxwellian, which 
is due to out-of-plane ion beam (localized current). 

We gather from panels (a)-(b) in Fig.([S|) (focus now 
only on solid curves) that in the strongly stressed X-point 
case, in about one Alfven time, super-thermal electrons 
show a full isotropy in all three spatial directions of the 
velocity distribution. In solar flare observations Kon- 



tar and Brown [l^ report that electron distributions are 
also nearly isotropic, which seems to contradict to what 
is expected from the standard model flare models that 
invoke strong downward beaming of electrons. In this re- 
spect, the match of our simulation results with the flare 
observations seem encouraging, in that stressed X-point 
collapse seems to be a viable mechanism acting during so- 
lar flares. For ions (panels (d)-(f) in Fig.®) behaviour is 
not so much different from the weakly stressed case (pan- 
els (d)-(f) in Fig.®), except for much higher velocities 
attained and distribution functions modified by kinetic, 
wave-particle interaction instabilities. The latter can be 
judged by sign changes in the slope of the distribution 
function, which can only occur when waves and particles 
exchange energy and momentum. 



IV. CONCLUSIONS 

By and large, the present work closes our initial study 
of stressed X-point collapse in the coUisionless regime 
started in Ref. [6| , by bridging gaps in the understanding 
of key physical aspects. The main findings can be listed 
as following: 

(i) despite significant differences of the initial setup 
between tearing unstable Harris current sheet [7[ and 
stressed X-point considered here, in both cases source 
of the reconnection out-of-plane electric field at the mag- 
netic null is provided by off-diagonal terms of the electron 
pressure tensor. 

(ii) we find that when mi/me S> 1 reconnection rate 
is independent of the ion-electron mass ratio and it is 
fast, which is also witnessed by Hesse et al. [§]• How- 
ever, when electron-ion mass ratio is unity, i.e. the Hall 
term is switched off, we show that reconnection rate is 
indeed slow. This broadly agrees with the results of Birn 
et al. (However see also Refs.[l3, [U [l3 for alterna- 
tive view). When the Hall physics is included, we also 
conjecture that the reconnection is fast because the mag- 
ntic field (being frozen into electron fluid, which moves 
significantly faster than ion fluid, as shown in Ref. [6|) 
is transported in and out of the diffusion region much 
faster than in the case of single fluid resistive MHD. We 



show that the amount of reconnected flux attained by 
LOcit — 25 in the cases of rrii/m^ ^ 1 and mi/nif. = 1 has 
the same ratio (w 4) as is the ratio of electron and ion 
speeds (w 4 — 5). 

(iii) we find that within one Alfven time, roughly ^ 
40% of the initial total energy (which is mostly stored in 
the magnetic field) is converted into the kinetic energy 
electrons, and somewhat more than half {^ 60%) into 
kinetic energy of ions. In solar flare observations a similar 
behaviour is seen |13j . 

(iv) When X-point is stressed strongly, in about one 
Alfven time, a full isotropy in all three spatial directions 
of the velocity distribution is seen for super-thermal elec- 
trons. Again similar behaviour is reported in solar flare 
observations 



14|. 



Resuming aforesaid, it seems that coUisionless, stressed 
X-point collapse is a viable mechanism for solar flares. 
Also, its behaviour is remarkably similar to tearing un- 
stable Harris current sheet which is thought to be more 
relevant for the Earth geomagnetic tail and generally to 
magnetospheric applications. 

We close this paper with some words of caution: 

(i) realising that size of the numerically simulated area 
40 X 40 electron skin depths (c/wpe = 0.16804 m in so- 
lar corona), i.e. 6.7 x 6.7 m is a tiny proportion of the 
actual solar flare site which can be tens of Mm. Thus, a 
priori it is not clear whether this mechanism is suitable 
to accelerate enough electrons in a much larger volume. 

(ii) whether e.g. the obtained electron-ion kinetic en- 
ergy partition (as the percentage of total released energy) 
being roughly 40% - 60%, would hold for a different mass 
ratio? We remind the reader that except in subsection 
III.B the ion-electron mass ratio was kept constant at 
mi/rrie = 100. We can only conjecture that this 40% 
- 60% partition holds for a different mass ratio, because 
we also found that the amount of reconnected flux versus 
time graph (Fig. 2) showed no mass ratio dependence. 
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